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AGILE  FREQUENCY  SELECTION  FOR  AN  ANTENNA  WITH  FREQUENCY 
DEPENDENT  BEAM  FIRING  ANGLE 


by 

H.C.  Chan 


Abstract 

(U)  This  work  addresses  the  potential  problems  and  their  solu¬ 
tion  of  employing  a  frequency  squinting  antenna  in  conjunc¬ 
tion  with  frequency  agility  in  a  naval  surveillance  radar 
system.  Frequency  agility  is  an  important  feature  in 
countering  enemy  jamming.  However,  when  it  is  applied 
with  a  frequency  squinting  antenna,  random  beam  firing 
angle  will  result.  This  random  beam  firing  angle  causes 
non-uniform  azimuthal  coverage.  Consequently  detection 
performance  in  certain  azimuths  scanned  by  the  antenna  will 
be  degraded.  An  optimization  technique  is  employed  to 
obtain  the  sequence  of  agile  frequency  selections  which 
provides  a  relatively  uniform  azimuthal  coverage,  while 
preserving  the  basic  random  properties  of  the  nominal 
sequence  of  agile  frequencies. 


1 .  INTRODUCTION 


(U)  In  modern  radar  systems  frequency  agility  is  an  important 

feature  to  counter  enemy  jamming.  The  radar  frequency  is  changed  from 
pulse-to-pulse,  or  from  burst-to-burst,  so  that  an  enemy  cannot  use  spot 
frequency  jamming  but  is  forced  to  use  less  effective  wideband  barrage 
noise  jamming  or  deploy  sophisticated  repeater  jammers  with  a  frequency 
following  capability.  With  pulse-to-pulse  frequency  agility,  the  selection 
of  the  instantaneous  radar  frequency  is  random  so  that  the  enemy  cannot 
predict  with  certainty  what  the  frequency  of  the  next  pulse  will  be.  The 
frequency  of  occurrence  of  the  radar  frequencies  is  roughly  uniform  so  that 
one  part  of  the  radar  band  is  not  favoured  over  another. 

(U)  With  a  conventional  (nonfrequency  squinting)  antenna  the 

selection  of  random  frequencies  is  straightforward,  since  the  look  direction 
is  independent  of  frequency.  This  is  not  the  case  for  frequency  squinting 
antennas,  however,  because  the  frequency  dependence  of  the  beam  direction 
can  cause  non-uniform  azimuthal  coverage  and  serious  degradation  in 
detection  performance. 
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(U)  Figure  1  depicts  the  geometry  of  the  frequency  squinting  antenna 

beam  firing  direction.  To  facilitate  the  analysis  of  this  problem,  the  fol¬ 
lowing  symbols  representing  the  relevant  parameters  are  defined: 

$  =  instantaneous  beam  firing  angle  with  respect  to  ship 
heading,  in  degrees 

0  =  bearing  of  the  antenna  normal  with  respect  to  ship 
heading,  in  degrees 

a  =  azimuthal  bias  angle  for  the  frequency  squinting 
antenna,  in  degrees 

Note:  This  is  the  beam  pointing  angle  with  respect  to  the 
antenna  aperture  normal  at  the  lower  limit  of  the  agile 
frequency  band.  For  travelling  wave  antennas,  boresight  beam 
condition  is  normally  avoided  because  of  the  existence  of  resonance. 

w  *=  antenna  rotational  rate  in  degree/sec 

ft  =  angular  velocity  of  the  antenna  normalized  to  the  radar 

pulse  repetition  interval  (PRI)  in  degrees/PRIfS  =  frequency 
squinting  sensitivity  in  degrees/GHz 
f  =  instantaneous  frequency,  in  GHz 

fj  *  lower  limit  of  the  agile  frequency  band,  in  GHz 
f2  =  upper  limit  of  the  agile  frequency  band,  in  GHz 
PRF  =  radar  pulse  repetition  frequency 

Y  =  3  dB  one=way  beamwidth,  in  degrees  of  the  frequency 
squinting  antenna  at  f^j 

(U)  A  comparative  simulation  of  the  azimuthal  coverage  between  a 

conventional  antenna  and  a  hypothetical  frequency  squinting  antenna  is 
performed.  Except  for  the  frequency  squint,  both  antennas  are  assumed 
to  have  identical  characteristics. 

(U)  At  time  t  =  0  the  normal  to  the  antenna  aperture  coincides  with  the 

ship  heading.  Let  T  be  the  time  interval  between  pulses  (PRI),  then 

T  =  1 /PRF  (1 ) 

(U)  A  pulse  will  be  emitted  whenever 

t  =  nT  n=0, 1 , 2, . . .  (2) 

(U)  The  index  n  will  be  called  the  pulse  index.  Since  the  angular 

velocity  of  the  antenna  is  normalized  to  the  system  PRI,  the  angular 
position  of  the  antenna  beam  at  the  time  of  pulse  emission  becomes  a 
function  of  the  pulse  index  n.  It  is  assumed  that  there  is  an  integer 
number  of  pulses  emitted  in  each  antenna  revolution.  The  beam  pointing 
direction  of  a  conventional  antenna  is  independent  of  radar  frequency,  hence 
it  is  governed  only  by  the  antenna  rotational  movement. 
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where  <pc { n )  denotes  the  angle  of  the  antenna  beam  with 
respect  to  the  ship  heading  at  the  time  the  nth  pulse  is 
emitted.  Subscript  c  denotes  a  conventional,  mechanically 
rotated  antenna. 

(U)  For  the  frequency  squinting  antenna,  the  beam  pointing  direction  is 
a  function  of  the  frequency  and  is  given  by: 


<t>s(n,f)  =  S(n)-[a  +  8(f  -  f  )] 

=  Qn-fa  +  6(f  -  f, )]  (4) 

where  a,  8,  £2  and  f^  are  defined  at  the  beginning  of  this 
section. 


4>3(n, f )  denotes  the  angle  between  the  antenna  beam 

and  the  ship's  heading  at  the  time  the  nth  pulse  of 
frequence  f  is  emitted.  Subscript  s  denotes  a 
mechanically  rotated,  frequency  squinting  antenna. 

The  beamwidth  of  both  antennas  will  be  affected  somewhat  by  the  radar 
frequency.  It  is  assumed  that  the  beamwidth  is  a  linear  function  of 
frequency  is  given  by; 


Y(f )  -  Y(f i )  -  5<f  -  f,)  (5) 

where  Y(f)  is  the  3  dB  one-way  beamwidth  at  frequency  f 
and  6  *  slope  of  the  change  of  beamwidth  with 
respect  to  a  change  of  frequency  (degree/GHz) 

This  has  been  observed  experimentally  [1]  to  be  a  good  approximation  for 
the  sandwich  wire  antenna  being  considered. 

(U)  For  the  purpose  of  the  simulation,  the  following  set  of  typical 

values  is  assigned  to  the  antenna  parameters 

Antenna  rotational  rate  oi  =  120°/sec 
PRF  =  740  Hz 

Frequency  squinting  sensitivity  8  =  24°/GHz 
ft  *  2.7  GHz 
f2  *  3.2  GHz 

Y  =  2.5°  @  2.7  GHz 
6  =  0.8°/GHz 


(U)  Comparisons  are  made  through  the  use  of  histograms  of  illumination. 

The  histogram  of  illumination  is  a  record  of  the  number  of  pulses  directed 
to  each  azimuth  in  one  antenna  revolution.  The  azimuth  is  divided  into  a 
number  of  equal  sectors.  For  each  pulse  transmitted,  one  or  more  of  the 
sectors  will  be  illuminated  by  the  3  dB  beamwidth.  Exactly  which  sectors 
are  illuminated  will  depend  on  the  beam  firing  direction  and  the  beamwidth 
at  the  chosen  frequency.  An  acumulator  is  assigned  to  each  of  the  azimuth 
sectors  which  is  initially  zeroed.  At  time  t  =  0,  the  antenna  normal  coin¬ 
cides  with  the  ship  heading.  The  first  pulse  is  emitted  at  this  time. 
Depending  on  the  3  dB  beamwidth  of  this  pulse  a  number  of  azimuthal  sectors 
in  the  vicinity  of  the  beam  direction  will  be  illuminated.  The  accumulator 
for  each  wholly  illuminated  sectors  will  be  incremented  by  one.  For  azimu¬ 
thal  sectors  which  are  only  partially  illuminated,  the  fraction  correspond¬ 
ing  to  the  illuminated  part  of  the  whole  sector  will  be  added  to  the  accumu¬ 
lator. 


THE  RANDOM  SEQUENCE 


(U)  The  agile  frequency  is  drawn  from  a  sequence  of  discrete  random  num¬ 
bers  which  assume  one  of  Nf  =  16  values  (0-15).  Each  value  will  corre¬ 
spond  to  one  of  16  frequencies  equally  spaced  across  the  agile  frequency 
band.  [The  restriction  of  the  frequency  division  to  discrete  values  of  Nf 
is  justified  because  the  radar  signal  has  a  finite  bandwidth).  The  separa¬ 
tion  between  the  frequencies  must  be  greater  them  the  signal  bandwidth  to  be 
significant.  If  the  number,  Nf  is  chosen  to  be  a  power  of  2,  one  can  uti¬ 
lize  the  so-called  pseudo-noise  (PN)  binary  sequence  [2].  For  example,  the 
grouping  of  continuous  binary  digits  in  twos  will  provide  4  values:  00=0, 
01=1,  10=2  and  11=3.  Grouping  in  fours  will  provide  16  values. 

(V)  Figure  2  shows  the  schematic  diagram  of  a  15-stage  PN  sequence 
generator.  This  PN  sequence  generator  will  produce  a  unique  binary 
sequence  of  length  2^-1  =  32767. 

(U)  Figure  3  shows  the  histogram  representing  the  relative  frequency 

of  occurrence  of  each  of  the  16  discrete  values  over  one  period  of  the  PN 
sequence.  Each  random  number  is  formed  by  four  contiguous  binary  digits. 
Figure  4  shows  the  auto-correlation  function  of  the  random  numbers.  The 
auto-correlation  function  is  defined  here  as: 


£  [x(n)-x] [x(n+k)-x] 
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Number  of  Occurrences  in  One  Period 


15-STAGE 


4  BITS  REPRESENTING 
A  UNIFORM  RANDOM 
NUMBER  WITH  16 
DISCRETE  VALUES 


MODULO  2  SUMMER 


FIGURE  2.  A  15-STAGE  PSEUDO-NOISE  GENERATOR. 


Frequency  Number 

FIGURE  3.  HISTOGRAM  OF  THE  UNIFORM  RANDOM  SEQUENCE. 


Tirp.fi  Lag  in  Number  of  Pulses 

FIGURE  4.  AUTOCORRELATION  FUNCTION  OF  UNIFORM  RANDOM  SEQUENCE 


where  x  =  mean  value  of  the  random  variable  x 


x 


1 

Nf 


i=0 


Nf  =  16 


(7) 


assuming  that  x(n)  are  uniformly  distributed  over  the  16  values  from  0  to 
15. 

(U)  It  can  be  seen  that  the  random  numbers  are  sufficiently  uniformly 

distributed,  and  they  are  essentially  uncorrelated  over  a  large  number 
of  samples.  This  ensures  that  it  would  be  difficult  for  an  enemy  to  predict 
the  frequency  of  the  next  pulse,  or  to  detect  any  bias  in  the  frequency 
selection. 


4.  AZIMUTHAL  COVERAGE  WITH  FREQUENCY  AGILITY 


(U)  For  a  conventional  reflector  type  antenna  operating  without  frequen¬ 

cy  agility,  the  number  of  pulses  on  target  per  scan  is  simply  given  by: 

03dBxPRF  2.3  x  740 

Hits  per  3  dB  beamwidth  =  -  =  -  =  14.183  (8) 

a)  1 20 

where  0  =  3  dB  beamwidth  of  antenna  at  mid  agile 

3dB 

frequency  band 

oj  =  antenna  rotation  rate  in  degrees/sec. 

(U)  Figure  5  shows  the  resulting  histogram  of  illumination  of  a  sector 
of  30°  over  one  antenna  revolution  by  the  squinting  antenna,  when  the  above 
random  sequence  is  used  to  select  the  agile  frequency.  The  term  histogram 
is  used  for  lack  of  a  more  descriptive  term.  It  should  be  emphasized  that 
it  is  not  a  probabilistic  parameter.  It  only  gives  a  measure  of  how  many 
pulses  are  to  be  received  in  any  particular  azimuth  over  one  revolution. 

This  histogram  is  typical  of  all  azimuthal  sectors.  It  shows  that  the  azi¬ 
muthal  coverage  is  no  longer  uniform  as  in  the  case  of  the  reflector  type 
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FIGURE  5.  TYPICAL  HISTOGRAM  OF  ILLUMINATION  OF  THE  SQUINTING 
ANTENNA  WITHOUT  AGILE  FREQUENCE  SELECTION 
OPTIMIZATION. 
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antenna.  (It  should  be  noted  that  the  use  of  this  PN  code  is  for  simulation 
purpose  only.  In  a  practical  system  it  would  be  prudent  to  employ  more 
secure  PN  codes  that  would  be  difficult  for  the  enemy  to  break  under 
tactical  conditions).  Certain  azimuths  receive  as  many  as  20  pulses,  while 
others  receive  as  few  as  8  pulses  over  one  revolution. 

( U )  For  adequate  detection  performance,  integration  of  pulses  is 

required.  In  the  case  of  pulse-to-pulse  frequency  agility,  signal 
integration  is  necessarily  noncoherent.  The  integration  gain  is  dependent 
on  the  number  of  pulses  available  for  integration.  Figure  6  [3]  shows  the 
required  signal-to-noise  ratio  (SNR)  for  0.95  probability  of  detection 
(PD)  as  a  function  of  number  of  pulses  integrated  by  a  linear  detector. 

These  curves 
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FIGURE  6.  REQUIRED  SNR  AS  A  FUNCTION  OF  THE  NUMBER  OF 
NON-COHERENT  PULSES  INTEGRATED.  FOR  A  LINEAR 
DETECTOR.  0.95  PROBABILITY  OF  DETECTION.  N0N- 
FLUCTUATING  TARGETS.  FOR  FIVE  VALUES  OF  Pf0 


are  for  nonfluctuating  targets,  and  for  five  values  of  false  alarm  probabil¬ 
ity  (Pfa)-  Assuming  for  the  moment  that  a  digital  sample  sorter  is  avail¬ 
able  to  rearrange  the  returned  signal  in  the  proper  azimuthal  order,  the  SNR 
for  each  azimuth  will  depend  on  the  number  of  pulses  received.  There  is  a 
difference  of  close  to  3  dBs  in  the  required  SNR  between  the  case  of  20 
pulses  integrated  and  the  case  of  8  pulses,  for  the  same  Pp  and  Pfa. 
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Taking  into  account  the  fact  that  there  is  a  difference  of  3  dBs  in  the 
energy  received  by  the  target  if  the  main  beam  is  exactly  on  target  compared 
to  the  case  when  the  edge  of  the  3  dB-beamwidth  frequeny  is  on  target,  there 
will  be  a  substantial  difference  in  the  detection  performance  from  one  azi¬ 
muth  to  another.  It  is  apparent  that  some  effort  must  be  devoted  to  finding 
the  proper  frequency  selection  algorithm  which  can  satisfy  the  frequency 
randomness  requirement  and,  at  the  same  time,  provides  relatively  uniform 
azimuthal  coverage. 

(U)  There  is  no  unique  solution  to  this  problem.  One  possible  approach 
is  to  use  optimization  techniques  which  is  the  topic  of  the  next  section. 


5.  CONSTRAINED  FREQUENCY  PERTURBATION  FOR  IMPROVED  AZIMUTHAL  COVERAGE 


(U)  In  order  to  obtain  a  uniform  coverage  when  the  squinting  antenna 

is  used  with  frequency  agility,  the  frequency  selection  must  be  modified. 
There  are  two  potentially  conflicting  requirements.  The  first  is  the 
requirement  of  uniformity  of  frequency  selection.  The  second  is  the  uni¬ 
formity  of  azimuthal  coverage.  A  compromise  must  be  struck  between  these 
two  requirements.  In  principle  one  can  start  out  with  a  sequence  meeting 
one  of  the  two  requirements  and  then  perturb  the  frequency  to  satisfy  the 
other.  It  is  highly  likely  that  an  attempt  to  improve  one  aspect  will  de¬ 
grade  the  other  but  by  observing  certain  constraints  in  the  frequency  per¬ 
turbation,  the  two  requirements  can  be  met  satisfactorily.  The  selection  of 
radar  frequency  affects  the  distribution  of  frequencies  as  well  as  their 
correlation.  These  two  measures,  however,  are  ensemble  averages;  that  is, 
the  distribution  and  correlation  of  frequencies  can  be  determined  only  after 
the  entire  sequence  is  determined.  Consequently,  it  is  more  logical  to 
start  with  a  sequence  of  frequency  selections  which  satisfy  the  distribution 
and  correlation  requirements  and  work  to  satisfy  the  uniform  azimuthal  cov¬ 
erage  requirement  by  frequency  perturbation. 

(U)  It  is  important  to  identify  the  objective  and  the  constraints. 

Since  the  antenna  rotational  rate  is  u>  and  the  pulse  repetition  frequency  is 
PRF,  the  total  number  of  pulses  to  be  emitted  per  antenna  revolution  is: 

Np  =  PRF  x  360/w  (9) 

(U)  If  the  azimuthal  coverage  is  divided  into  Np  sectors,  then  the 

angular  extent  of  each  sector  will  be: 

A0  =  360°/Np  (10) 

(U)  The  most  uniform  coverage  would  result  if  one  and  only  one  pulse  is 

directed  to  the  centre  of  each  of  the  Np  azimuthal  sectors. 


Consequently  the  objective  is  to  manipulate  the  frequency  in  order  to 
achieve  the  above  condition.  This  may  be  considered  as  a  constrained 
optimization  problem.  One  of  the  constraints  to  be  observed  is  the 
resulting  distribution  of  the  frequencies  selected.  It  must  be  reasonably 
uniform.  If,  during  the  course  of  modifying  the  frequency  selection, 
certain  frequencies  are  selected  significantly  more  often  than  others,  it 
would  result  in  a  bias  in  the  frequency  distribution.  This  could  be 
exploited  by  enemy  jammers.  Another  constraint  is  the  correlation  of  the 
frequencies.  There  appears  to  be  no  direct  control  which  can  be  applied  to 
guarantee  an  adequately  uncorrelated  sequence.  The  undesirable  effects  can 
be  minimized  however  by  keeping  the  number  of  changes  made  to  the  original 
frequency  selection  to  a  minimum. 

(U)  Figure  7  is  a  tabulation  of  the  number  of  pulses  pointed  to 
identical  azimuthal  sectors  in  one  antenna  revolution.  Each  sector  is  A0  in 
angular  extent.  For  this  particular  random  sequence,  the  frequencies 
selected  resulted  in  792  azimuthal  sectors  to  which  no  pulse  is  directed. 


N  =  Number  of  Pulses  Received 
FIGURE  7.  DISTRIBUTION  OF  MULTIPLE  PULSE  CELLS. 


There  are  821  azimuthal  sectors  which  have  exactly  one  pulse  directed  to 
them.  There  are  457,  119,  27  and  4  azimuthal  sectors  which  received  two, 
three,  four  and  five  pulses  respectively  during  one  antenna  revolution.  In 
subsequent  discussion,  in  order  to  simplify  the  description,  the  term  "empty 
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cell"  will  be  used  to  describe  the  azimuthal  sectors  which  did  not  receive 
any  pulse  in  one  antenna  revolution.  By  "receiving",  it  is  meant  that  the 
peak  of  the  antenna  beam  falls  into  the  azimuthal  extent  of  that  particular 
sector.  The  term  “multiple  pulse  cell"  will  be  used  to  describe  the 
azimuthal  sectors  which  receive  more  than  one  pulse  during  the  course  of  one 
antenna  revolution.  The  extra  pulses  of  the  multiple  pulse  cells  can  be 
diverted  to  neighbouring  azimuthal  sectors  by  a  slight  increment  or 
decrement  of  the  radar  frequency.  This  can  be  accomplished  without 
disturbing,  to  any  significant  degree,  the  inherent  randomness  of  the 
original  sequence.  The  redirection  of  a  pulse  to  other  azimuthal  sectors  by 
changing  its  frequency  will  be  called  "beam  switching". 

(0)  The  large  fluctuation  of  azimuthal  coverage  is  due  to  congregation 
of  empty  cells  and  multiple  pulse  cells.  Table  1  tabulates  the  number  of 
consecutive  empty  cells  in  one  antenna  revolution.  If  cells  nearby  have 


(U)  Table  1  Distribution  of  consecutive  empty  cells 


No.  of  consecutive  empty  cells 
revolution 

No.  of  occurrence  per  antenna 

1 

315 

2 

115 

3 

51 

4 

13 

5 

3 

6 

2 

7 

0 

8 

1 

multiple  pulses,  these  empty  cells  can  be  filled  by  diverting  one  of  the 
extra  pulses  from  the  multiple  pulse  cells  with  a  slight  change  in  its 
frequency.  However,  there  is  a  moderate  number  of  empty  cells  which  are 
joined  together.  For  example,  there  is  a  string  of  eight  azimuthal  sectors 
which  are  not  illuminated.  Since  the  ability  of  the  antenna  to  back  scan  is 
limited  by  the  squinting  range,  it  is  highly  likely  that  a  number  of  the 
empty  cells  cannot  be  covered  by  directly  switching  the  extra  pulses  in  the 
multiple  pulse  cells. 

(U)  Assuming  that  the  antenna  is  rotating  in  the  counter-clockwise 
direction  and  that  the  ship  heading  is  designated  as  the  true  north  in  the 
radar  display,  the  beam  firing  angle  will  always  be  lagging  behind  the 
antenna  normal.  The  converse  is  true  if  the  antenna  is  rotating  in 
the  clockwise  direction  or  the  signal  is  being  fed  from  the  opposite  end  of 
the  antenna.  For  the  purpose  of  illustration,  the  convention  depicted  in 
Figure  1  will  be  used. 


S.S.'f;- ,*■  •»  r%  -  v  v  y\-  y  vVv  y ^  f  j '  v  • '» - 


£ 


i 


(U)  Suppose  in  one  antenna  revolution,  the  jth  azimuthal  cell  has 

received  N  pulses  with  frequencies  fj,  f2»  f3»...,  fN,  respectively, 

then  N-1  of  these  pulses  can  be  diverted  so  that  they  can  be  pointed  to  em¬ 
pty  cells.  Suppose  the  (j-l)th  cell  is  an  empty  cell,  in  order  to  direct 
one  of  the  extra  pulses  to  it,  the  frequency  of  that  pulse  must  be  changed 
by  the  amount 


Af  =  A0 

T 


where  3  is  the  squinting  sensitivity  and 
A0  is  defined  in  (10) 

(U)  If  the  angular  distance  between  the  multiple  pulse  cell  and  the 

empty  cell  is  too  great,  then  the  frequency  increment  required  may  cause  the 
agile  frequency  range  to  be  exceeded.  This  difficulty  can  be  alleviated  by 
involving  all  the  pulses  directed  to  azimuthal  cells  located  between  the  em¬ 
pty  cell  and  the  multple  pulse  cell.  For  instance,  let  the  jth  cell  be  a 
multiple  pulse  cell  and  the  (j-k)th  cell  be  an  empty  cell.  There  are  k  azi¬ 
muthal  cells  separating  the  empty  cell  and  the  multiple  pulse  cell.  The 
change  in  frequency  required  to  direct  an  extra  pulse  from  the  jth  cell  to 
the  (j-k)th  cell  is  therefore; 


Af*  =  kAf  =  kA0 


Suppose  now  that  the  frequency  of  the  pulse  in  question  can  only  be  incre¬ 
mented  by  Af  without  exceeding  the  agile  frequency  range.  Consequently,  it 
cannot  be  used  to  directly  fill  the  empty  cell;  however,  it  can  be  used  to 
point  the  beam  to  its  immediate  neighbor,  the  (j-1 )th  cell,  using  a  frequen¬ 
cy  increment  of  Af.  The  pulse  originally  assigned  to  the  (j-l)th  cell  can 
then  be  directed  to  the  (j-2)th  cell  by  introducing  the  same  amount  of  fre¬ 
quency  increment.  This  process  continues  until  the  empty  cell  is  filled. 


(U)  It  is  possible  that  even  this  procedure  will  fail  to  fill  all  the 

empty  cells  because  of  the  frequency  agility  band  limitation.  Furthermore, 
this  procedure,  if  applied  indiscriminately,  will  destroy  the  uniformity 
property  of  the  frequency  distribution.  Frequencies  at  the  agile  frequency 
band  edges  would  tend  to  occur  more  often.  Thus  a  better  procedure  is  need¬ 
ed  to  complete  the  task. 

(»J)  In  the  following,  a  set  of  steps  is  outlined  which  will  result  in  a 

uniform  histogram  of  illumination  for  one  antenna  revolution  while,  at  the 
same  time,  attempting  to  maintain  random  frequency  assignments. 
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A  nominal  frequency  selection  sequence  is  determined  from  a 
discrete  uniform  random  number  generator.  The  selection  has 
Np  entries.  The  frequencies  are  indexed  and  stored  in  an 
array  so  that  the  kth  entry  represents  the  frequency  used 
for  the  kth  pulse. 

A  coarse  histogram  of  illumination  is  constructed  based  on 
the  nominal  sequence  by  computing  the  exact  direction  of  the 
beams.  The  circumference  of  the  radar  coverage  is  divided 
into  Np  sectors,  each  with  an  azimuthal  width  given  by 
(10).  The  cells  are  indexed  in  ascending  order  with  the 
first  cell  representing  the  azimuthal  cell  in  the  ship’s 
heading.  There  is  an  accumulator  assigned  to  each  of  the 
azimuthal  cell.  If  the  beam  pointing  direction  of  a  pulse 
falls  in  one  of  the  cells,  the  accumulator  assigned  to  this 
cell  is  incremented  by  one. 

The  frequencies  of  the  pulses  are  then  adjusted  slightly  so 
that  the  beam  will  be  directed  precisely  at  the  centre  of 
the  cell. 

A  beam  destination  table  is  also  constructed  by  recording 
the  index  of  the  azimuthal  cell  to  which  a  pulse  is  direct- 


The  coarse  histogram  of  illumination  is  scanned  for  multiple 
pulse  cells.  Once  a  multiple  pulse  cell,  say  the  kth  cell, 
is  located,  the  cells  immediately  before  and  after  the  kth 
cell  will  be  examined  to  determine  if  it  is  an  empty  cell. 

If  it  is,  then  one  of  the  pulses  originally  assigned  to  the 
kth  cell  will  be  switched  to  it  by  changing  its  frequency  by 
an  amount  of  ±  Af  given  by  (11).  The  plus  sign  is  taken 
when  the  empty  cell  is  lagging  the  kth  cell,  while  the  nega¬ 
tive  sign  is  taken  if  the  empty  cell  is  leading  the  kth  cell 
in  azimuth.  The  resulting  frequency  is  checked  to  ensure 
that  it  does  not  exceed  the  agile  frequency  band  limits, 
the  steps  are  repeated  until  all  empty  cells  have  only  sin¬ 
gle  pulse  neighbours.  The  beam  destination  table  and  the 
histogram  of  illumination  are  updated  after  each  beam 
switching  has  taken  place. 

The  resulting  coarse  histogram  of  illumination  is  again 
scanned  for  multiple  pulse  cells.  When  one  is  located,  the 
closest  empty  cell  is  found.  It  is  then  filled  by  switching 
the  pulse  originally  directed  to  its  immediate  neighbouring 
cell  located  between  the  empty  cell  and  the  multiple  pulse 
cell.  Thus  the  empty  cell  effectively  propagates  towards 
the  multiple  pulse  cell.  This  continues  until  the  empty 
cell  appears  immediately  next  to  the  multiple  pulse  cell 
which  can  then  be  filled  by  assigning  one  of  the  extra 
pulses  to  it.  These  steps  are  repeated  until  all  empty 
cells  are  filled  except  those  which  cannot  be 
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reached  without  employing  frequencies  exceeding  the  agile 
frequency  band  limits. 

(U)  Steps  a  to  f  will  be  called  preliminary  beam  switching  and  they 

do  not  involve  any  optimization. 

g.  After  the  above  steps  are  performed  the  histogram  of  illumi¬ 
nation  will  be  fairly  uniform  except  for  a  small  number  of 
spots  where  empty  cells  cannot  be  eliminated  by  simple  beam 
switching.  This  may  be  due  to  the  large  angular  distance 
between  the  multiple  pulse  cell  and  the  empty  cell,  or  it 
may  be  due  to  the  fact  that  the  frequencies  of  some  pulses 
are  very  close  to  the  agile  frequency  band  edge  and  cannot 
be  incremented  or  decremented  any  further.  In  any  case,  the 
solution  would  involve  switching  beams  selectively. 


(U)  At  this  point  it  would  be  useful  to  introduce  an  objective  function 
which  is  a  function  of  the  frequency  histogram.  The  frequency  histogram  is 
a  tally  of  the  number  of  occurrences  of  each  of  the  Nf  discrete  frequen¬ 
cies  over  one  antenna  revolution.  Because  of  the  way  the  discrete  random 
numbers  are  obtained,  the  nominal  frequency  histogram  is  fairly  uniform  ex¬ 
cept  for  the  smallest  frequency  value.  After  preliminary  beam  switching, 
the  uniformity  of  the  frequency  distribution  will  be  degraded.  The  ensuing 
steps  are  designed  to  determine  the  proper  switching  of  pulses  so  that  the 
remaining  empty  cells  can  be  eliminated  while,  as  far  as  possible,  the  uni¬ 
formity  of  the  frequency  histogram. 


<U) 


The  objective  function  is  defined  as: 


1  vf  2 

-  >  (NF  -  F  ) 

Nf  i=l  i  M 


(13) 


where  NF^  is  the  number  of  occurences  of  the  ith 


discrete  frequency,  and 
Fm  is  the  mean  value  of  the  frequency  of 
occurrence 


Fm  =  Np/Nf 


(U)  The  procedures  outlined  below  may  be  considered  locally  optimal  in 
the  sense  that  they  are  concerned  with  finding  the  best  way  to  switch  a 
pulse  originally  directed  to  a  particular  multiple  pulse  cell  to  a  particu¬ 
lar  empty  cell.  Consider  a  multiple  pulse  cell  i  and  an  empty  cell  l. 
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Assume  that  there  are  K  azimuthal  cells  separating  them  as  depicted  in  Fig¬ 
ure  8.  Depending  on  the  number  of  intermediate  cells  to  be  involved  in  the 


Azimuthal 


K  Azimuthal  Calls 
Separating  call  fi  and  calif  l 


FIGURE  8.  MULTIPLE-STAGE  BEAM  SWITCHING  STRATEGY. 


beam  switching,  there  are  a  large  number  of  alternative  ways  in  which  the 
empty  cell  can  be  filled.  For  instance,  it  may  be  possible  to  first  divert 
one  of  the  N  pulses  which  are  originally  directed  to  the  multiple  cell  i  to 
the  jth  cell,  and  then  change  the  frequency  of  the  pulse  for  the  jth  cell  so 
that  it  will  point  to  another  intermediate  cell  k.  Finally  the  frequency  of 
the  pulse  for  cell  k  is  changed  so  that  it  will  now  point  to  the  empy  cell 
A.  The  criterion  for  selecting  a  particular  way  of  beam  switching  is  that 
it  minimizes  the  objective  function  (12).  This  problem  can  be  solved  using 
dynamic  programming  [4]  for  multiple  stage  decision  processes.  This 
approach,  however,  requires  a  large  amount  of  memory  because  it  is  a  delayed 
decision  process.  All  relevant  items  must  be  temporarily  stored  as 
alternative  decisions  are  attempted  at  each  stage.  Alternatives  which 
result  in  a  higher  objective  function  are  purged,  thereby  reducing  the  total 
memory  requirement.  It  is  much  more  effective  than  direct  enumeration  which 
attempts  all  possible  alternatives  and  then  selects  the  optimum.  For  the 
problem  at  hand,  the  dynamic  programming  approach  becomes  very  cumbersome 


when  the  number  of  intermediate  beam  switches  is  large.  If,  on  the  other 
hand,  the  number  of  intermediate  beam  switches  is  limited  to  one,  then  the 
problem  becomes  a  2-stage  decision  process.  In  this  case  the  dynamic  pro¬ 
gramming  approach  is  feasible  and  can  be  handled  with  direct  enumeration. 
Aside  from  practical  considerations,  there  are  other  reasons  to  limit  the 
problem  to  a  two-stage  beam  switching  problem.  The  first  is  the  fact  that 
the  majority  of  empty  cells  can  be  reached  by  pulses  originally  assigned  to 
a  multiple  pulse  cell  via  only  one  intermediate  cell.  The  second  is  that  in 
order  to  minimize  the  detrimental  effects  of  beam  switching  on  the  random 
properties,  the  number  of  frequency  changes  should  be  kept  to  a  minimum. 
Finally  in  the  rare  occasions  that  an  empty  cell  cannot  be  filled  with  only 
one  two-stage  beam  switch,  then  two  or  more  two-stage  beam  switches  can  be 
cascaded  to  complete  the  task. 

(U)  An  example  should  aid  to  illustrate  the  operations  of  these  proce¬ 
dures.  Assume  that  azimuthal  cell  i  is  a  multiple  pulse  cell  and  cell  l  is 
an  empty  cell.  The  problem  is  to  find  a  way  of  switching  one  of  the  extra 
pulses  originally  directed  to  cell  i  to  cell  l.  There  are  many  alternative 
ways  in  which  this  can  be  accomplished.  Figure  9  symbolically  shows  just 
two  of  such  alternatives. 

(U)  For  alternative  No.  1,  the  beam  switching  sequence  is  as  follows: 

i.  The  pulse  originally  directed  to  cell  i  is  redirected  to  cell 

j* 

ii.  The  pulse  originally  directed  to  cell  j  is  redirected  to  cell 

l. 

(U)  For  alternative  No.  2  the  beam  switching  sequence  is  as  follows: 

i.  The  pulse  originally  directed  to  cell  i  is  redirected  to  cell 

k. 

ii.  The  pulse  originally  directed  to  cell  k  is  redirected  to  cell 

l. 

(U)  The  selection  is  based  on  the  inflicted  cost  between  choosing  cell  j 
or  cell  k  as  the  intermediate  cell. 

(U)  In  addition  the  following  items  are  given: 

a.  the  frequency  of  the  extra  pulse  originally  pointed  to  cell 
i  is  f  ^ . 

b.  The  frequency  of  the  only  pulse  pointed  to  cell  j  is  f2. 

c.  If  the  pulse  originally  pointed  at  cell  i  is  to  be  switched 

to  cell  j,  its  frequency  must  be  changed  to  f3. 

d.  If  the  pulse  originally  pointed  to  cell  j  is  to  be  switched 
to  cell  l,  its  frequency  must  be  changed  to  f4* 

e.  The  frequency  of  the  only  pulse  pointed  at  cell  k  is  fg. 

f.  If  the  pulse  originally  pointed  at  cell  is  to  be  switched 

to  cell  k,  its  frequency  must  be  changed  to  fg. 

g.  If  the  pulse  originally  pointed  at  cell  k  is  to  be  switched 
to  cell  l,  its  frequency  must  be  changed  to  £7. 


Beam  Switching  Alternative  No.  1 


Pulse  Originally  Pointed  To 
Cell  #1  Is  Redirected  To 
Cell  ff]  By  Changing  Its 
^  Frequency  To  fj 

Cellfi  \ 


(Multiple  Pulse  Cell) 


Pulse  Originally  Pointed  To 
Cell  ft  Is  Redirected  To 
Cell  jfk  By  Changing  Its 
Frequency  To  fg 


Pulse  Originally  Pointed  To 
Cell  f\  Is  Redirected  To 
Cell  fl  By  Changing  its 
Frequency  To 

N Cell  § l 
(Empty  Cell) 


Pulse  Originally  Pointed  To 
Cell  #k  Is  Redirected  To 
Cell  f  By  Changing  Its 
Frequency  To  fj 


Beam  Switching  Alternative  No. 2 


The  Costs  For  The  Two  Alternatives  Are 
Computed,  And  The  One  With  The 
Smaller  Cost  Is  Selected 


FIGURE  9.  TWO-STAGE  BEAM  SWITCHING  STRATEGY. 
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(U)  If  the  first  alternative  is  chosen,  the  count  of  the  frequency 

histogram  for  fi  will  be  reduced  by  one,  and  the  count  for  f3  will  be 
increased  by  one.  Similarly,  the  count  for  f 2  will  be  decreased  by  one, 
and  the  count  for  f4  will  be  increased  by  one.  The  resulting  objective 
function  due  to  this  sequence  of  beam  switching  is: 


_f  2 

)  (NF  -  F  )  +  (NF  -  1 

i*1  1  M  1 

i*2  +  (NF  +  1 

3 

i*3  +  (NF  -  1 

2 

i*4  +  (NF  +  1 

4 


(U)  If  the  alternative  No.  2  is  chosen,  the  count  for  f.j  will  be 

reduced  by  one  and  the  count  for  fg  will  be  increased  by  one.  Simi¬ 
larly  the  count  for  fg  will  be  decreased  by  one  and  the  count  for  f 7 
will  be  increased  by  one.  The  resulting  objective  function  due  to  this 
sequence  of  beam  switching  is 

,.f  2  2 

F  =  )  (NF  -  F  )  +  (NF  -  1  -  F  ) 

2  i*1  i  M  1  M 

2 

i*6  +  (NF  +  1  -  F  ) 

6  M 

2  (15) 

i*5  +  (NF  -  1  -  F  ) 

5  M 

2 

i*7  +  (NF  +  1  -  F  ) 

7  M 
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The  choice  between  these  two  alternatives  is  determined  by  selecting  the  one 
which  yields  the  smaller  objective  function  F.  The  optimum  policy  is  found 
by  considering  all  possible  intermediate  cells  such  as  cell  j  and  k  in  the 
above  example.  The  selected  alternative  will  be  replaced  whenever  another 
alternative  yields  a  smaller  objective  function  F.  These  optimization 
procedures  continue  until  all  empty  cells  in  the  coarse  histogram  of 
illumination  are  filled. 


KSS 


(U)  A  set  of  subroutines  is  written  which  implement  these  procedures 
together  with  a  program  which  demonstrates  their  use. 


(U)  Figure  10  shows  the  detailed  histogram  of  illumination  resulting 
from  a  frequency  selection  obtained  by  this  program.  Trace  A  is  the 
coverage  before  optimization.  Trace  B  shows  the  azimuthal  coverage  after 
optimization.  It  can  be  seen  that  the  resulting  histogram  of  illumination 
is  very  uniform.  The  total  number  of  pulses  received  at  any  azimuth  varies 
within  plus  and  minus  one  pulse. 
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FIGURE  10.  COMPARISON  OF  AZIMUTHAL  COVERAGE  OF  AGILE 
FREQUENCY  SELECTIONS  WITH  AND  WITHOUT 
OPTIMIZATION. 


(U)  It  remains  to  check  the  random  properties  of  the  optimized  frequency 

selection.  After  optimization  the  frequencies  obtained  are  quantized  to  the 
nearest  of  the  16  values  available.  Figure  11  is  the  frequency  histogram  of 
the  optimized  sequence.  As  can  be  seen  it  is  fairly  uniform.  Figure  12 
shows  the  autocorrelation  function  of  the  optimized  frequency  selection  com¬ 
puted  from  the  quantized  frequencies.  It  shows  no  increase  in  correlation 
among  frequencies. 

(U)  The  optimized  results  for  azimuthal  coverage  obtained  above  (Figure 

10)  assume  that  the  radar  frequency  can  be  changed  continuously.  It  would 
be  of  interest  to  find  what  the  degradation  would  be  in  the  azimuthal  cover¬ 
age  if  the  radar  frequency  is  quantized.  Figure  13,  Figure  14  and  Figure  15 
show  the  resulting  histogram  of  illumination  when  the  optimized  frequency 
selection  is  quantized  to  16,32  and  64  discrete  levels  respectively.  It  can 
be  seen  that  at  64  levels,  the  azimuthal  coverage  is  quite  acceptable. 
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FIGURE  11.  COMPARISON  OF  FREQUENCY  HISTOGRAMS  OF  THE 
NOMINAL  AND  THE  OPTIMIZED  AGILE  FREQUENCY 
SELECTIONS. 
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FIGURE  12.  AUTOCORRELATION  FUNCTION  OF  THE  OPTIMIZED 
AGILE  FREQUENCY  SELECTION. 
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FIGURE  15.  TYPICAL  HISTOGRAM  OF  ILLUMINATION  FOR  THE 
OPTIMIZED  AGILE  FREQUENCY  SELECTION 
QUANTIZED  TO  64  DISCRETE  LEVELS. 


6.  PERIODICITY 


(U)  The  sequence  of  agile  frequency  selection  obtained  with  the  above 

procedures  contains  Np  entries,  which  is  the  total  number  of  radar  pulses 
emitted  per  antenna  revolution.  Thus  by  the  time  the  antenna  completes  one 
revolution,  the  sequence  would  have  been  exhausted.  The  same  sequence  can 
be  repeated  and  the  uniformity  in  azimuthal  coverage  will  be  assured. 
However,  this  creates  a  periodicity  in  the  frequency  selection.  If  an  enemy 
ESM  receiver  possesses  some  memory  capability,  it  is  conceivable  that  this 
periodicity  can  be  exploited.  This  is  not  an  insurmountable  problem.  One 
way  to  reduce  the  periodicity  is  to  have  several  random  sequences  each  of 
which  will  provide  uniform  azimuthal  coverage.  These  sequences  will  be 
constrained  such  that  they  are  identical  in  their  entries  at  the  beginning 
of  the  sequences.  This  will  enable  the  sequences  to  be  switched  from  one  to 
the  other  without  losing  the  uniformity  of  azimuthal  coverage.  The 
switching  of  the  sequences  can  be  done  randomly,  thereby  making  the  apparent 
period  extremely  long.  Nevertheless,  this  problem  must  be  considered  when 
evaluating  the  performance  of  a  practical  system. 
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7.  CONCLUSION 


(U)  The  use  of  frequency  agility  in  a  radar  system  with  a  mechanically 

rotated,  frequency  squinting,  antenna  causes  a  random  beam  pointing  direc¬ 
tion  dependent  on  the  radar  frequency  employed.  This  random  beam  pointing 
direction  interacts  with  the  mechanical  rotational  movement  of  the  antenna 
and  causes  a  non-uniform  distribution  of  the  transmitted  pulses  over  all 
azimuths.  As  a  result,  detection  performance  may  differ  significantly  from 
one  azimuth  to  another.  For  a  simple  naval  surveillance  radar  system,  this 
problem  can  be  alleviated  by  employing  optimization  techniques  in  the  selec 
tion  of  the  sequence  of  agile  frequencies.  Simulation  results  showed  that, 
with  the  agile  frequency  selected  from  32  or  64  discrete  frequencies  spaced 
equally  across  the  agile  frequency  band,  the  azimuthal  coverage  can  be 
brought  to  within  ±2  and  ±1  pulses,  respectively,  of  that  for  the  conven¬ 
tional  antenna  where  the  average  number  of  pulses  per  beam  is  14.  The  dif¬ 
ference  in  detection  performance,  with  one  pulse  more  or  less  integrated, 
would  have  a  neglible  effect  on  performance.  The  same  technique  can  be 
employed  to  obtain  the  agile  frequency  selection  in  the  burst-to-burst  mode 
when  MTI  filtering  is  used. 

(U)  Although  the  frequency  squint  does  not  seem  to  present  an  insur¬ 

mountable  problem,  it  does  represent  an  added  constraint  on  the  signal  pro¬ 
cessing  system.  Consequently,  the  effect  of  the  frequency  squint  on  the 
overall  signal  processing  system  performance  must  be  analyzed  when  a 
specific  system  is  implemented. 
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